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Abstract 

An angular analysis for B decays into two vector mesons gives parameters sensitive to new 
physics. They are the decay amplitudes of the three helicity states, the asymmetry in triple product, 
and the CP violating phase. The measurements of these quantities are performed for B° — > J/ipK*° 
decays in a data sample with 253 fb _1 taken by the Belle detector in the KEKB asymmetric 
e + e~ collider. With the time-integrated angular analysis, the decay amplitudes' moduli and phase 
angles are measured to be \A \ 2 = 0.585 ± 0.012 ± 0.09, |Aj|| 2 = 0.233 ± 0.013 ± 0.008, \A ± \ 2 = 
0.181 ±0.012 ±0.008, arg(Aj|) = 2.888 ±0.090 ±0.008 radians, and arg(A_L) = 0.239 ±0.064 ±0.010 
radians. The difference between arg^y) and arg(^4j_) is 2.649 ± 0.110 radians, which is shifted 
from 7r by more than 4cr; this can be interpreted as evidence for the existence of the final state 
interaction. The amplitude values are used in triple product correlations to obtain asymmetries 
of A { f = 0.101 ± 0.033 ± 0.007 and A { ft = -0.102 ± 0.032 ± 0.003 for B° meson decays, and 
Aft = 0.056 ± 0.030 ± 0.006 and Aft = -0.091 ± 0.028 ± 0.003 for W meson decays. The 
time dependent angular analysis gives the measurements of CP violating parameters sin20i and 
cos 20i to be 0.30 ± 0.32 ± 0.02 and -0.31 ± 0.91 ± 0.10, respectively. The value of cos20i is 
—0.56 ± 0.86 ± 0.11 if sin 2^1 is fixed at the world average value (0.731). 

PACS numbers: 13.25.Hq, 14.40.Nd 
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INTRODUCTION 



An angular analysis of B meson decay to two vector mesons can be a sensitive probe 
of new physics. Since the analysis is equivalent to the simultaneous study of three decay 
modes corresponding to the possible helicity states of two vector mesons, the effect of the 
new physics that might appear in the interference among the states can be observed cleanly 
with the cancellation of most systematic effects. 

Recently, the anomaly in the decay amplitudes in B — > <j)K* has been reported where 
the longitudinal amplitude is lower than the Standard Model prediction[l|. For the precise 
study of such anomalies, a comparison with theoretically clean reference modes is necessary. 
The decay B° — > J/ipK*° is an ideal mode for this purpose since it is a tree dominated mode 
with a very low penguin contribution where the theoretical ambiguity on the CP-sensitive 
parameter sin 20x0] is less than 1% level[3]. 

There are three classes of parameters that can be used to probe new physics obtained 
through the angular analysis. The first is the measurement of the decay amplitudes of the 
three helicity states. They can be obtained by the time-integrated angular analysis to flavor 
specific decays. The comparison of the amplitudes between flavor-separated samples probes 
the direct CP violation. The second is the triple product correlation, which can be extracted 
from the measured decay amplitudes. This quantity is sensitive to the T-violation. The third 
class is comprised of the CP parameters (sin20i and cos20i) that are measured through the 
time-dependent angular analysis. In particular, the measurement of cos 20! , which appears 
in the time-dependent interference terms, is important both to solve the two-fold ambiguity 
in 20x and to test the consistency of this determination with the more precise value from 
other b — > ccs decays. 

In this paper, we report the measurements of all three classes of parameters for B° — > 
J/ipK*° decays. 



DATA SAMPLE 

The data sample used in this analysis corresponds to an integrated luminosity of 253 fb _1 
recorded with the Belle detector at the KEKB electron-positron collider Events are 
required to satisfy the hadronic event selection criteria |2| and have R2 < 0.5, where R2 is 
the ratio of the second to zeroth Fox- Wolfram event shape moments 0. 

The reconstruction of J/ip candidates is performed using the dilepton decays J/ip — > 
e + e~ and fi + fi~. For the e + e~ mode, a track is identified as an electron or positron by 
a comparison of the energy measured in the electromagnetic calorimeter (ECL) with the 
momentum measured in the central drift chamber (CDC), the shape of the cluster energy 
deposit in the ECL, the specific ionization (dE/dx) measured in the CDC, and the light 
yield in the aerogel Cerenkov counters (ACC). The likelihoods that the track is an electron 
or a hadron are determined from these measurements. The likelihood ratio is required 
to be consistent with the electron hypothesis. To correct for energy lost by final state 
radiation, the energy of any cluster in the ECL within 0.05 radians of the initial track 
momentum is added to that of the track. The invariant mass of each so-augmented electron- 
positron pair is calculated, and the pair is associated with J/ip — > e + e~ if the mass is 
in the range 2.95 GeV/c 2 < M(e + e~) < 3.15 GeV/c 2 . Tracks are identified as muons 
rather than hadrons by means of likelihoods based on (i) a comparison of the number of 
layers with associated hits in the muon detector (KLM) with the number expected based 
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on momentum and (ii) the energy of the associated hit in the ECL. An oppositely charged 
pair of muons is identified as arising from J/ip — > fi + fi~ if the invariant mass is in the 
range 3.05 GeV/c 2 < M(/i + /i~) < 3.15 GeV/c 2 . To improve the momentum resolution, a 
kinematic fit that uses the J/ip mass as a constraint is performed on J/ip candidates passing 
the above selections. 

Candidate K*° mesons are reconstructed using the decay modes K*° — > K + ir~ and 
K*° -> K%7T°. A track is identified kaon based on the ratio of its likelihoods to be a 
kaon vs a pion. These likelihoods are obtained from the measurements mentioned above from 
the TOF, CDC and ACC detectors. Tracks that are not identified as kaons and not used 
as leptons in the J/ ip reconstruction are treated as charged pion candidates. Kg candidates 
are reconstructed from pairs of oppositely charged tracks that satisfy three conditions: 1) 
the distance of closest approach of each track to the nominal interaction point is larger than 
0.03 cm, 2) the angle between the Kg momentum vector and the vector displacement of 
the Kg vertex point from the J/ip vertex is less than 0.15 radians, and 3) the reconstructed 
decay vertex of the Kg is at least 0.1 cm away from the interaction point. Each pair with 
an invariant mass satisfying 0.47 GeV/c 2 < M(7r + 7r~) < 0.52 GeV/c 2 is identified as a Kg 
meson. To improve the mass resolution, the pion momenta are re-fitted, constraining both 
tracks to originate at the reconstructed Kg vertex. 

Candidate 7r° mesons are reconstructed from clusters in the ECL that are unmatched to 
charged tracks and having energy greater than 40 MeV. A photon pair with an invariant 
mass in the range 0.125 GeV/c 2 < M(77) < 0.145 GeV/c 2 is identified as a tt°. A mass- 
constrained fit is performed to obtain the momentum of the 7r° meson. A Kit pair is 
considered a K* candidate if the invariant mass of the pair is within 74 MeV/c 2 of the 
nominal A'*(892) mass. 

Candidate B° mesons are reconstructed by selecting events with a J/ip and a K* meson 
and examining two quantities in the center-of-mass of the T(4S): the beam-constrained 
mass (Mbc) and the energy difference between the B candidate and the beam energy (AE). 
The beam-constrained mass, which is the invariant mass of a reconstructed J/tp and K* 
calculated taking the energy to be the beam energy, is required to be in the range 5.27 — 
5.29 GeV/c 2 . For the mode with K*° — > K + ir— , \AE\ is required to be less than 30 MeV, 
while AE is required to satisfy -50 MeV < AE < 30 MeV for K*° -> K^n°. 

To eliminate slow tt° backgrounds, the angle 9 k* of the kaon with respect to the K* 
direction in the K* rest frame is required to satisfy cos9k* < 0.8. This is equivalent to 
demanding that the ir° momentum be greater than 175 MeV/c. When an event contains 
more than one candidate passing the above requirements, the combination for which both 
Mbc and AE are the closest to the B° mass and zero, respectively, is selected. After all the 
selections, 8194 K*° — > K + 7t~ candidates and 354 K*° — > A^vr candidates remain. Figure 
^ shows the distributions of M^ c and AE for the samples. 



MEASUREMENT OF DECAY AMPLITUDES 

The decay amplitudes of B — > J/ipK* decays are measured in the transversity basis 0]. 
The definition of the angles is shown in Fig. El The direction of motion of the J / ip in the 
rest frame of the B candidate is defined to be the x-axis. The y-axis is chosen along the 
direction of the projection of the K momentum into the plane perpendicular to the x-axis 
in the B rest frame. The x-y plane contains the momenta of the J/ip, the K, and the 7r 
mesons. The z-axis is then perpendicular to the x-y plane according to the right-hand rule. 
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FIG. 1: AE and Mb c distributions for reconstructed B° — > J/ipK*°(K + 7T~) (upper figures) and 
£?° — > J/-0i^*°(X^7r°) (lower figures) samples. Solid lines show the results of two dimensional fits 
while dashed lines are estimated background contaminations. 

The angle between the the positive lepton daughter and the z-axis in the J/ip rest frame 
is defined as tr . The angle between the x-axis and the projection of the l + momentum onto 
the x-y plane is defined as <pt r in the same frame. The angle 9k* is defined as described in 
the previous section. 

The distribution of these three angles for B — > J/ipK* decays is described in terms of 
three amplitudes ||: 

1 dT 9 

o~l a — 7T~ = * cos 2 9 K *(1 - sin 2 9 tr cos 2 (j) tr )\A \ 2 

1 a cos 9 tr a cos 9 K *d<pt r o2tx 

+ sin 2 9 K * (1 - sin 2 9 tr sin 2 tr ) \A\\ | 2 

+ sin 2 9 K * sin 2 9 tr \A ± \ 2 

+r] sin 2 9k* sin 29 tr sin tr Im(A*\A±) 

j= sin 29k* sin 2 9 tr sin 2(f) tr Re{AQA\\ ) 

+rj—= svo.29k* sm.29 tr cos0t r Jra(-A o \Aj_)], (1) 
v2 

where Aq, A\\ and v4j_ are the complex amplitudes of the three helicity states in the trasver- 
sity basis, and 77 = +1 (—1) for B° (B°). The coefficient |v4 | 2 denotes the longitudinal 
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I- 

FIG. 2: The definition of transversity angles. 

polarization of J/ip while \A±\ 2 (|^4||| 2 ) gives the transverse polarization component along 
the z-axis (y-axis). In the decay mode where K*° — > K^tt , \Aq\ 2 + |^4||| 2 (|^4_i_| 2 ) is the 
component corresponding to the CP-even (CP-odd) state. 

The amplitudes are determined by fitting this function to the measured three-dimensional 
distribution in 9 tr , (p tr and 9 k*, taking into account the detection efficiency and background. 
The resolution of the angular measurements is estimated by Monte Carlo (MC) simulation 
and found to be typically less than 0.02 radians. The value of rj is determined from the 
charge of pions or kaons used in the K* reconstruction. We use only B° — > J/i/jK*°(K + tt~) 
decays for the measurement since r\ is not well-defined in B° — ► J /ipK* (K^ii®) decays. 

The fit is performed using an unbinned maximum likelihood method. The probability 
density function (Q) is defined using the theoretical distribution in (JTJ) and can be expressed 
as 

1 d 3 T 

G(x,y,z,M hc ,AE) = N x [f sl (M hc ,AE) x e(x,y,z) x - (x,y,z) 

I axayaz 

+ Y^f l c f (M hc ,AE)xA cf (x,y,z) 

i 

+ f nr {M hc ,AE) x A nr\X, y, Z) 

+ f cb {M hc , AE) x Acb{x, y, z)}, (2) 

where x = cos9 tr , y = <ptr, and z = cos9k*, N is the normalization factor of Q, e(x,y,z) is 
the detection efficiency as a function of the three angles, and A c f, A nr , A cb are the angular 
shapes for the cross-feed, non-resonant and combinatorial backgrounds, respectively. 

Here, f sig (M hc , AE), p cf (M bc ,AE), f nr (M hc ,AE), and f cb (M hc ,AE) are the fractions 
of signal events, cross-feed contamination, non-resonant contamination and combinatorial 
background, respectively, as a function of Mt, c and AE. These fractions are obtained with 
a data subsample of 140 ftT 1 by the procedure given in Ref. with an improvement to use 
both Mb c and AE in the fits. 

The detection efficiency function e(x, y, z) is obtained from a large MC sample of 6 million 
events generated polarized with the decay amplitudes set at the values of the previous 
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TABLE I: Summary of measured decay amplitudes. 





Flavor averaged 


B° BO 


IA,| 2 


0.585 ±0.012 


0.581 ±0.016 0.589 ±0.016 


IV 


0.233 ± 0.013 


0.220 ±0.018 0.246 ±0.018 


\Ai_\ 2 


0.181 ±0.012 


0.199 ±0.018 0.164 ±0.017 


arg(A||) 


2.888 ± 0.090 


2.937 ±0.136 2.854 ±0.120 


arg(A ± ) 


0.239 ± 0.064 


0.303 ±0.090 0.182 ±0.089 



measurement [9] . Events are histogrammed in a 20 x 20 x 20 grid in the cos6 tr — 4>tr — cosOk* 
cube. The distribution is fitted to the product of three one-dimensional polynomials. 

The angular distribution function for the cross-feed background (A c f) is determined from 
MC simulation. The function for the non-resonant production (A nr ) is determined from 
events in the sideband of the Kn mass distribution (1.0GeV/c 2 < M(Kn) < 1.3GeV/c 2 ). 
The distribution for the combinatorial background (A c b) is obtained from events in the side- 
band with 5.2 GeV/c 2 < Mb c < 5.26 GeV/c 2 . These normalized distributions are parameter- 
ized as the product of three one-dimensional polynomials whose parameters are determined 
from the fit. 

In the fit, the phase of is defined to be zero relative to those of the other amplitudes 
since the overall phase of the decay amplitudes is arbitrary. The other five parameters, 
\Aq\ 2 , |v4||| 2 , |^4±| 2 , axg(^4||) and axg(A±) are free parameters in the fit. The normalization 
condition of the amplitudes 

|A | 2 ±|A||| 2 ±|AJ 2 = 1 (3) 
is taken into account by adopting the extended likelihood defined as 

N obs 

- In L = - £ In Gi + N exp - N obs \n(N exp ) (4) 

8=1 

where N b s is the number of events used for the fit. N exp is defined to be N b s ■ ( | ^4 1 2 + |^.]| | 2 + 
|t4jJ 2 ) to incorporate the normalization condition. The normalization N of Q is recalculated 
by numerical integration whenever the fit parameters change. The parameter values deter- 
mined from the fit are summarized in Table and the projected angular distributions are 
shown in Fig. El The distributions are corrected for the effects of detector acceptance and 
backgrounds. Small discrepancies from or ir are observed in arg(A||) and &rg(A±). The 
difference of these two phases is 2.649 ±0.110 radians which is shifted from 7r by more than 
4<7. This is read to be evidence for the existence of final state interaction. 

There are two choices for the imaginary phases of the helicity amplitudes 0). We obtain 
the above values by taking the conventional choice of the phases noted as solution I in 
Ref.[l(j. If the alternative set of the phases is used, the imaginary phases become arg(A||) — > 
- arg(A||), and arg(A ± ) ->• tt - arg(A ± ). 

The parameter values are also measured for B° (4121 candidates) and B° (4073 candi- 
dates) samples separately and the results are given in Table The measured values are 
consistent between both B° and B° decays, indicating no evidence for direct CP violation. 

Systematic uncertainties in the fit are determined as follows: 1) efficiency function (MC 
statistics and effect of polarization), 2) angular distribution functions for backgrounds, 3) 
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FIG. 3: Distributions of projected angles. Solid lines show results of the fit. The data points are 
corrected for the detector efficiency and the backgrounds are subtracted. 



TABLE II: Systematic errors in the measurement of decay amplitudes. 



Item 


|A)| 2 |^||| 2 |^4±| 2 arg(A||) arg(4i) 


Efficiency 
PDF for backgrounds 
Background fractions 
Slow pion efficiency 
Polarization in NR 


0.003 0.002 0.001 0.001 0.003 
0.002 0.002 0.002 0.004 0.006 
0.001 0.002 0.003 0.001 0.002 
0.008 0.007 0.007 0.006 0.007 
0.002 0.002 0.003 0.003 0.003 


Total 


0.009 0.008 0.008 0.008 0.010 



background fractions, 4) slow pion efficiency, and 5) polarization effect in non-resonant 
decays. These contributions to the systematic error are summarized in Table |TI] The 
dominant contribution arises from the uncertainty in the detection efficiency for slow pions. 
This is estimated by comparing the results with a cut at cos 9k* < 0.9. The angular shape 
in non-resonant backgrounds is a potential source of uncertainty since it might contain 
tails of other resonances whose polarization is different. The uncertainty is determined by 
comparing with the results obtained assuming the angular distribution to be that of phase 
space decay. 
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TRIPLE PRODUCT CORRELATIONS 

The T- violating triple product correlations are sensitive probes to new physics effects 
The asymmetry of the triple product is defined as 

= £(g • (g x g) > 0) - T{v[ ■ (g x g) < 0) 

T r(ui • (?r 2 x ^ 3 ) > 0) + r(«i • x u 3 ) < 0) u 

where ul is the momentum of or K*° meson, while v<i and Vs are the polarization vectors 
of J/ijj and meson, respectively. Experimentally, two asymmetries can be defined using 
the measured decay amplitudes for B° decays, as follows: 



Im(A ± A*) (2) _ HUg 



The corresponding asymmetries for B° decays are defined as and A^ . The Standard 
Model predicts tiny values for these asymmetries and no difference between B° and B° 
mesons. 

By putting the measured amplitude values in Eq. |H1 the triple product asymmetries are 
obtained as 

A? = 0.101 ±0.033 ±0.007 
A (2) = -0.102 ±0.032 ±0.003 
A? = 0.056 ±0.030 ±0.006 
A {2) = -0.091 ±0.028 ±0.003 

As seen, all the obtained triple product asymmetries are close to zero, furthermore, no 
difference between A T 's and A T 's is observed. 



MEASUREMENT OF CP VIOLATION PARAMETERS 

The CP violation parameters sin 2<fti and cos 2<p 1 are measured by studying the angular 
distributions in the CP decay mode B — > J/ipK*°; K*° — > K^tt as a function of the decay 
time difference between B° and B° (At). The At-dependent form of Eq. ^becomes0: 

d 4 r(9 tr ,cf> tr ,9 K *,At) e -l^l/- s o e 

2^gi{0tr,(Ptr,9K*)a>i(At) (7) 



d cos 9t r d(ptrd cos 6k*(1 At 2tb» 

where tb is the lifetime of a B° meson, and the angular terms are expressed as: 

gi = 2 cos 2 9k* (1 — sin 2 9 tr cos 2 4> tr ) (8) 
g 2 = sin 2 9 K *(l - sin 2 9 tr sin 2 <j) tr ) (9) 
g 3 = sin 2 9k* sin 2 9 tr (10) 

g>4 = —= sin 29k* sin 2 9 tr sin 2<p tr (11) 
V2 

#5 = sin 2 9k* sin 29 tr sin (p tr (12) 

g»6 = —= sin 26*^. sin 29 tr cos tr (13) 
V2 
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while the amplitude terms are expressed as: 

a x = \A \ 2 (1 + sin 20i sin AmAt) (14) 

a 2 = |v4||| 2 (l ± sin 20 x sin AmAt) (15) 

a 3 = \A ± \ 2 (1 - r/ sin 20i sin AmAt) (16) 

a 4 = Re(A* { A )(l ± r? sin 20 x sin AmAt) (17) 

a 5 = T)Im(AT,A±) cos Am At - r]Re(A*A ± ) cos20! sin Am At (18) 

a 6 = ?7/m(AoA ± ) cos AmAt - r]Re(A* A ± ) cos20! sin AmAt. (19) 

Here, Am is a B° — B° mixing parameter, rj is +1 for B° and —1 for B°. A®, An and A± are 
the decay amplitudes given in the previous section. Two CP violation parameters appear in 
the formula, viz, sin 20! and cos 20! . The latter parameter appears in the interference terms 
a 5 and a%. 

The procedures to measure the decay time difference and to determine the flavor of decay- 
ing B° meson are described elsewhere jy]. The values of sin20i and cos20i are determined 
by fitting the function to the measured angles and At, taking into account the detection 
efficiency and background. The fit is done using an unbinned maximum likelihood method. 
The probability density function for an event is defined as 

V = f si g{M hc ,AE)e{6 tr ,<j) tr ,6 K »)- 

a cos Vtrd<Ptrd cos Ok* d/Xt 

-|A*|/t b o 

+^ {E f l cf(M bc , AE)A cf (6 tr , <p tr , 9 K .) 

+ fnr(M hc , AE)Anr(9tr, <j>tr, 9 K *)} 

+6(At)f cb (M bc , AE)A cb {e tr , cf> tr , 6 K *). (20) 

Here, f S i g , f c f, f nr and f c b are the fractions of signal, cross feed, non-resonant production, 
and combinatorial background components, respectively, while A c f, A nr and A c b are the 
corresponding three-dimensional angular shape functions. Also, e is a three dimensional 
detection efficiency function for the signal. These functions are the same as those used in 
the decay amplitude measurement. The flavor tagging procedure gives the flavor q of tag- 
side B meson and the probability w that this flavor determination is incorrect. To account 
for the effect of wrong tagging in this fit, 77 is replaced with —q(l — 2w) within the signal 
angular distribution in Eq. |2H1 

Each term in V is then convolved with the appropriate resolution functions separately 
for the signal, backgrounds having the B° lifetime (namely, cross feeds and non-resonant 
production) and the combinatorial background with a zero-lifetime 5-function shape, which 
takes into account the effect of vertex position smearing caused by the detector resolution, 
charmed meson decays, and poorly reconstructed tracks. The resolution parameters are 
calculated event by event using the energy of the reconstructed B meson, the direction of 
the B meson with respect to the beam axis and the momentum of the B meson in the 
center-of-mass frame 2 ■ 



In the fit, the helicity amplitudes are fixed at the values determined in the previous 
section. The values of the lifetime and mixing parameter are set at the values given in 
Ref. ^3[. The CP parameters sin20i and cos20i are the only free parameters in the fit. 
From the fit to the data, we obtain 

sin 20i = 0.30 ±0.32 ±0.02, 
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TABLE III: Systematic errors in the measured CP parameters. 



Item 


sin 2cj)i 


cos 2(pi 


cos 20i (sin 2</>i = 0.731) 


Vertexing 


0.008 


0.050 


0.051 


Resolution parameters 


0.005 


0.010 


0.010 


Wrong tagging fractions 


0.009 


0.061 


0.062 


^-function for BG 


0.006 


0.004 


0.005 


tb and Am 


0.004 


0.010 


0.010 


Decay amplitudes 


0.005 


0.041 


0.048 


Angular shape for backgrounds 


0.007 


0.025 


0.044 


Background fraction 


0.011 


0.021 


0.021 


Total 


0.020 


0.096 


0.107 



cos 20i = -0.31 ±0.90 ±0.10. 

When we fix the value of sin 20i at the world average value (0.731)jl3, the value of cos 20i 
becomes 

cos20! = -0.56 ±0.86 ±0.11 (21) 

We obtain the above values by taking a set of the phases given in the previous section. If 
the alternative set of the phases is used in the fit, the sign of cos20i flips, giving cos20i = 
+0.31 ± 0.86 ± 0.11, while sin20i does not change. 

Systematic uncertainties in the fit are determined in the same manner as those in b — > ccs 
sin 20i measurement 0] . In addition, the uncertainties that come from the angular analysis 
(decay amplitudes, background angular shapes, and background fractions) are estimated. 
The uncertainty in decay amplitudes are determined by varying the values by one standard 
deviation. The uncertainty in the background fractions and angular shapes are estimated 
in the same manner as that in the decay amplitude measurement. 

The raw asymmetry in the measured At between samples tagged as B° and B° is shown 
in Fig. |3] together with the result of the fit. 

The systematics in the fit are checked by applying the same fitting procedure to the 
sample of B —>■ J/tpK*°(K + 7r~) decays. We obtain following values for parameters sin20i 
and cos 20i : 

"sin 20i " = -0.037 ±0.068 
"cos 20i " = 0.103 ±0.179 

These are consistent with zero, as expected. 
CONCLUSION 

A full angular analysis is performed for B° — > J/i/jK*° (K*° — > K + n~) decays. The 
complex decay amplitudes of three helicity final states are measured by the simultaneous 
fit to three transversity angles and the best fit values are \A \ 2 = 0.585 ± 0.012 ± 0.009, 
|A||| 2 = 0.233±0.013±0.008, \A ± \ 2 = 0.181±0.012±0.008, arg(An) = 2.888 ±0.090 ±0.008 
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-5 5 
At(ps) 

FIG. 4: Raw asymmetry in measured At between samples tagged as B° and B°. The solid line 
shows the projection of the fit. 

radians, and aig(A±) = 0.239 ± 0.064 ± 0.010 radians, where arg(A ) is defined to be zero. 
There are two choices for the set of imaginary phases. We obtain the above values by taking 
the conventional choice of phases. The measured value of \Aj_\ 2 shows that the CP even 
component dominates in B° — > J / ip K*° ; K*° — > Kgir decays. Small discrepancies from 7r 
or are observed in arg(A||) and arg(Aj_), respectively. The difference of these two phases 
is calculated to be 2.649 ± 0.110 radians and the shift from n is more than 4o\ This is 
read to be an evidence for the existence of final state interactions. The amplitudes are also 
measured for B° and B° decays separately. The best fit values are consistent between them 
and no direct CP violating effect is observed. 

The measured amplitudes are converted into the asymmetries in the triple product corre- 
lations: the best fit values are Aft = 0.101 ±0.033±0.007 and A% ] = -0.102±0.032±0.003 
for B° mesons, and A$p = 0.056 ± 0.030 ± 0.006 and A^ ] = -0.091 ± 0.028 ± 0.003 for W 
mesons. The differences between two asymmetries for B° and B° mesons are consistent with 
zero and no new physics effects are observed. 

The time dependent angular analysis is performed for B° — > J/ipK*° (K*° — > Kgir ) 
decays. By the simultaneous fit to the three transversity angles and At, the CP violation 
parameters sin 20! and cos 20! are determined to be sin 2<pi = 0.30±0.32±0.02 and cos 20! = 
— 0.31±0.91±0.10. Fixing sin 20! at the world average value (0.73) gives cos 20! = — 0.56± 
0.86 ± 0.11. The decay amplitudes are fixed at the values measured in the time integrated 
fit with the conventional choice of imaginary phases. When taking the other choice, the sign 
of cos 20i flips; however, we cannot put a constraint on the sign of cos 20! regardless of the 
choice of the imaginary phase set. 
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TABLE IV: Comparison of obtained decay amplitudes and CP violation parameters with previous 
measurements. Here, the sign of cos 20i of Belle measurements is put to be consistent with BaBar's 
choice of phases. 





This Measurement 


Bellef9]ri6] 


BaBar[15] 


Luminosity 


253fb _1 


29fb 1 (78 for CP) 


83fb (113 for CP) 


\M 2 

\Ai_\ 2 
arg(A||) 
&ig{A L ) 


0.585 ±0.012 ±0.009 
0.233 ±0.013 ±0.008 
0.181 ±0.012 ±0.008 
2.888 ± 0.090 ± 0.008 
0.239 ±0.064 ±0.010 


0.62 ±0.02 ±0.03 

0.19 ±0.02 ±0.03 
2.83 ±0.19 ±0.08 
-0.09 ±0.13 ±0.06 


0.566 ±0.012 ±0.005 
0.204 ±0.015 ±0.005 
0.230 ± 0.015 ± 0.004 
2.729 ±0.101 ±0.052 
0.184 ±0.070 ±0.046 


sin 20i 
cos 201 


0.30 ±0.32 ±0.02 
±0.31 ±0.91 ±0.11 


0.13 ±0.51 ±0.06 
±1.40 ± 1.28 ±0.19 


-0.10 ±0.57 
+3.32lg;^ ± 0.27 



Table IIVI shows the comparison with previous measurements. The measured decay 
amplitudes are consistent among all measurements. 
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